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SUMMARY 

A curved member under axial load is analyzed using beam column theory to 
determine nonlinear response and the tangent stiffness associated with small 
displacements from the nonlinear state. Such a result is suitable for incor- 
poration into a general nonlinear analysis using a corotational coordinate 
system to describe the rigid body type motion of individual members. Applica- 
tion of the method to buckling problems is also shown. Several examples are 
given to show the accuracy of the method. 

INTRODUCTION 

The nonlinear analysis of frame or lattice structures can be accomplished 
with 3 variety of approaches with varying degrees of success depending on ihe 
problem, in the present paper* attention is given to frame structures that may 
undergo large displacements but small ;rain which is often characterized as a 
very flexible or lightly loaded structure. The approoch taken is to use a 
moving or corotational coordinate system that follows each member so there is 
no limit on displacement and rotation. Exact solutions of the beam column 
equations are used for total and incremental forces so that accuracy is 
achieved with few or no member interior nodes. Such a theory has been 
developed sn reference 1 for straight members. For very flexible structures, 
member imperfections may be a significant factor in nonlinear response or a 
member may he <w*signed to have a small curvature. In the present paper, the 
method of reference 1 is extended to include curved members. The application 



of the tangent stiffness matrix associated with the incremental forces to 
buckling problems is also illustrated. Results of nonlinear analysis of a 
simple arch and imperfect columns are compared with previous work to indicate 
the accuracy of the approach. 


SYMBOLS 


Ai,A2»A3,A4 

B 


C 

bi»b2,b3,b4 

ci,C2 

Cb 


0,d 


dp 

EA 

El 

e 

F.f 

6 


Gi,G2 

H 

K.k 


constants of integration in equation (A4) 
transformation matrix between relative displacements (or 
associated forces) and incremental displacements (or associated 
forces) in local coordinates 
matrix defined by equation (A21) 

bowing functions defined by equations (A9) through (A12) 
stability functions defined by equations (A6) 
shortening due to bending 

global and local displacement vectors respectively 

local relative displacement vector 

axial stiffness of member 

bending stiffness of member 

amplitude of member imperfection 

global and local force vectors respectively 

force vector associated with relative displacements 

matrix giving contribution of total member loads to incremental 

local force vector 

bending extension coupling terms defined by equation (A18) 
extensional stiffness term defined by equation (A18) 
global and local stiffness matrices respectively 
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L member length 

L* member length after deformation 

mi,m2 bending moments at ends of member 

p external load on structure In X direction 

PE euler load of column 

Pl»P2 axial loads at ends of member 

Q ^ (mj + m2)/L* 

q ratio of axial load to Euler load 

qx>q 2 shear loads at ends of member 

R,S stiffness matrices from reference 5 

T transformation matrix from local coordinates of displaced 

member to global coordinates 

U. u global and local displacements In X(x) direction 

V, v global and local displacements In Y(y) direction 

W component of T matrix 

w deflection of member 

X,Y global coordinate system 

x,y local coordinate system 

B Imperfection parameter defined by equation (AS) 

Y angle between X axis and line passing through member ends 

0 angle between axial and diagonal member 

01, 02 relative rotations at ends of member 

tp rotation at a node 

4> stability function parameters defined In equation (A 6 ) 

A stability function parameters defined In equation (A 6 ) 



Subscripts 


a 

d 


i.j 

0 


denotes axial member 
denotes diagonal member 
denote nodes 1 and j 

Initial value corresponding to zero displacement 


Superscripts 

T denotes matrix transpose 

' denotes derivative with respect to q 

A 6 preceding a variable Indicates Incremental value 


ANALYSIS APPROACH 

The equations developed are for plane frames but extension to space frames 
as was done In reference 2 Is straightforward. Consider a member in a plane 
frame connecting node 1 to j. After displacement, Its position Is Indicated 
in figure 1. The analysis problem. Illustrated In figure 2, is to determine 
the forces pj, mj, m 2 , when the member has undergone the relative deformation 
of 01 , 02 and U 2 « Note that p 2 , qi, d2 are determined by equilibrium 

from pi, mi, m 2 . In addition, the tangent stiffness relating Incremental dis- 
placements to Incremental forces must be determined. 

6f « kfid (1) 

where 


6fT 


®di » 

6mj , 5 p 2 * 

4 ^ 2 * 

6m2] 

6dT 

« [ 6ui, 

6vi , 

6(>ii, 6 u2, 

6V2 , 

4 <I' 2 ] 


4 



ORIGINAL PAGE 13 
OF POOR QUALITY 


The matrix k can be assembled Into a global stiffness by the transformation 


where 


K * T^kT 


T » 


W 0 
0 W 


( 2 ) 

(3) 


W 


cos r sin r O' 

-sin Y cos Y 0 

0 0 1 


where y is defined by 


cos Y = (Xj - Xi + Uj - Ui)/L* 
sin Y = (Yj - + Vj. V,)/L* 

with 

L* = v'lXj - XT .Uj - Ui)^ V(Yj . . Vy- V^)-2 


(4) 


(5) 


In assembling the global stiffness matrix, the forces in the members can be 
combined with any external load to give the unbalanced force at each node. The 
negatives of these forces are applied to the structure to achieve a new deflected 
position and the process is repeated until convergence is achieved. At each 
stage of the iteration, the relative displacements are determined as shown in 
figure 1 by 


* i'l • (Yi - Yo) 
62 = - (Yj - Yo) 

U2 = L* - L 
Ui = 0 
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where Yq and L ar the value of y and L* respectively evaluated for 
* Uj » Vi = Vj » 0, The matrix k and expreslons for mj, m2 PI 
derived In reference 1 for straight members. The corresponding development 
for curved members Is given In Appendix A. 

If a member load Is high relative to Its buckling load and/or the member 
has undergone large relative deformations (Oj, 62 large), the approach 
described above may become Inaccurate. One remedy is to put Internal nodes in 
the member. However, this can greatly increase solution times as a result of 
the increases in matrix size but even more significant Is the increase in band 
width. An alternative approach is to use substructures for individual members, 
as indicated in figure 3. An arbitrary number of internal nodes is located in 
the member; three internal nodes are shown in figure 3 for an example. The 
analysis of the single member is performed by the method described above for 
the displacement boundary conditions shown in the figure. This will require 
several iterations to obtain the correct set of internal displacements. If 
the nodes are numbered starting in the interior so that the ends are numbered 
last, the tangent stiffness matrix for the converged displacements may be tri- 
angulated down to the last two nodes. The remaining matrix is the tangent 
stiffness matrix for the entire member. The forces at the ends of the member 
can be determined to complete all the information required for the next 
iteration in the analysis of the complete structure having nodes only at con- 
nections between members. 
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RESULTS 


The present theory has been applied to several problems to assess the 
accuracy of the analysis. In the following sections, the elastica problem, 
column with initial imperfection, the deformation of a shallow arch, and buckling 
of a lattice truss with member imperfections will be presented. 

The Elastica 

The elastica problem has a known solution valid for unlimited deflection 
and so provides a good basis for comparison. Results presented in Table 1 were 
calculated for the geometry of figure 4 and are compared with the exact solu- 
tion which is in the form o^ elliptic Integrals from reference 3. One half of 
a pin ended column was modelled with a symmetry boundary condition in the 
middle. Results are presented for no internal nodes and for one internal node 
by use of the substructure analysis. The comparison with the exact result is 
good to load levels appreciably higher than the buckling load, with the results 
for one internal node essentially identical to the exact result up to twice the 
buckling load. Note the large rotations present that correspond to relative 
rotations well over one radian. 

Imperfect Column 

The same model as in figure 4 was analyzed with an imperfection equal 
.OIL. The results are shown in figure 5 where the end load is plotted as a 
function of the lateral deflection. For loads less than the buckling load, 
the result is in close agreement with the classical amplification factor formula. 
Above 30% of the buckling load, the two results deviate with the present solution 
approaching the elastica for higher loadings. The use of a substructure with 
one internal node shows a slight difference from the result without a 
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substructure at loads above the buckling load Indicating that reasonable 
accuracy has been achieved without an Internal node. 

Shallow Arch 

A clamped shallow arch with a central load has been analyzed and compared 
with finite element results froiTi reference 4, In figure 6, the load P Is 
plotted as a function of the central deflection U. One node at the point of 
load application was used. Results without substructuring and with a sub- 
structure having one Internal node are In close agreement Indicating convergence 
has been achieved. This Is confirmed by the close agreement with the finite 
element results of reference 4. Note that 32 curved elements with three 
degrees of freedom per node were required for convergence. At the higher 
loadings where the arch has become Inverted, the present results indicate the 
finite element solution may not be completely converged. 

Buckling of a lattice Truss with Imperfect Diagonal Members 

In reference 5 an analysis was made of tne buckling of a family of lattice 
configurations with repetitive geometry. A feature of the method Is that the 
six degrees of freedom of one node are all that is Involved In the solution. 
Figure 7 Is an Illustration of two of the configurations analyzed. The present 
theory allows the treatment of configurations with imperfect members. The 
details of how the analysis of reference 5 Is extended to Include Imperfect 
members Is given In Appendix 6. 

An example of the application of the theory Is given In figure 8. A 
three element truss, three bays long and simply supported at the ends has been 
analyzed. Note that a truss with a length that Is a multiple of three bays 
would also buckle at the same load unless the column mode was lower. The 
buckling load for a perfect truss Is presented as a function of length In 
8 


i 



reference 5. The diagonal members Mere assumed to have an Imperfection as 
shown on the figure. The buckling load, normalized by the load that would be 
achieved If the diagonals provided simple support Is plotted as a function of 
the amplitude of the Imperfection. The symbols are from a complete analysis uf 
the structure by the computer program of reference 6 which Is a general space 
frame buckling and vibration program based on exact member theory for straight 
members. The Imperfection was modelled by two straight segments with a maxi- 
mum deviation from a straight line equal to the amplitude of the Imperfection. 
Such a model Is somewhat stiffen than the curved shape used In the present 
analysis which accounts for the slightly higher result. However, the close 
agreement throughout the large range of Imperfection Indicates the accuracy 
of the present approach. For this problem, the Imperfection has a very large 
effect on the buckling load. 


CONCLUDING REMARKS 

The tangent stiffness matrix for a curved member based on the exact folu- 
tlon of the beam column equation has been developed. The matrix would be suit- 
able for incorporation Into a nonlinear analysis that uses a corotatlonal 
coordinate system. The accuracy of the method as shown by several numerical 
examples Is such that no Internal nodes would be required in any member for 
many applications. In cases where the relative rotation Is large or the axial 
load Is approaching or exceeding the Individual member buckling load a substruc- 
turing method Is presented which allows the complete analysis to be accomplished 
with nodes only at points of load application or at connections to other members. 
The only limit to such an approach Is the excedence of the elastic limit. 

The tangent stiffness matrix can also be used for calculating buckling loads 
for certain configurations having Imperfect members. 
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Applications of the method are shoim for several problems to ladicate the 
accuracy of the approach. The problems treated are (1) the elastics, ( 2 ) ar 
Imperfect column, (3) centrally loaded shallow arch, and (4) buckling of a 
lattice truss with Imperfect members. The benefit of using exact member equa- 
tions compared to the usual finite element approximation functions was demon- 
strated clearly with the arch problem. The present method with only one 
node at the point of load application achieved results comparable to a solution 
using 32 curved finite elements In conjunction with a corotatlonal coordinate 
system. 
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APPENDIX A 

NONLINEAR STIFFNESS OF A CURVED MEMBER 

The expressions for member stiffness are based on an exact solution of the 
beam column equation for a member with a sinusoidal deviation from a straight 
line. The axial force In the member Is calculated accounting for shortening 
due to member curvature. The differential equation for lateral deflection Is 

EIw*'^ + Pi(Wo“ + w") « 0 (Al) 

with boundary conditions 

w(0) ■ w(L) * 0 

w'(0) « ej w'(L) « 62 (A2) 

and 

Wq * sin wx/L (A3) 

1 he solution for w Is 

w - Aj A 2 x/L + A 3 cos ^x/L + A 4 sin <>x/L + 3^s1n irx/L (A4) 

1 -q 

where q Is the ratio of the applied compressive load to the Euler load 

q = PiL2/(»2ei) 
ai.J 

4 “ 

The bending moments at each end are found by the usual process to be 
mi » E1 [ci(6i-6) C2(62'*^S)3A 

(A5) 

m 2 ■ EI[c2(6i-0) + ci(02+6)]A 
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where 


*■” A 

and c^, C 2 are the stability functions defined as 

cj ■ A (sin ♦ - cos i) C 2 • A (1-cos 

♦ 


(A6) 


where 


A • 2(1 - cos ♦) - s jn ^ 

♦ 


Note that cos ^ and sin ♦ are real for t Imaginary so the above expressions are 

♦ 

valid for negative q. The compressive force p^ Is given by 


Pi • -EA( U 2 ■ ui ♦ C5) 


(A7) 


where 




bl(®l'*’®2^^ + b2(9i-62“2$)2 + b 2 ( 0 j- 62 ) + b^0 


(A8) 


and 


00 ■ e»/L 
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The bj and 63 are bowing functions derived b> Onan in reference 2 who 
showed the following relations with the stability functions 


b 

b 


1 

2 







(A9) 

(AlO) 


where a prime denotes derivative with respect to q. 

The b 3 and b^ are new bowing functions due to the Initial deflection and 
<ire given by 


b 

3 




1 ] 


Peflne a loading vector fp as 


fp^ • [Pi* <"1. in2] 


and a displacement vector dp as 


dp^ ** lu2 - uj, «j, Ojj] 


Incremental values of fp and dp may be related by 


6fp • t iSdp 

The components of the matrix t are obtained as 


(All) 

(A12) 


(A13) 


(A14) 


(Alb) 


t 


U 



r 1 

r.l 


♦ iVl J3. 

3q adpj 


(Alb) 
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(A17) 


where 

H = 1 + EAL^[bj(6j + 62)^ “ ®2 * ^ ^3(®1 “ ®2^®o * 

61 * -[2bi(ei + 02) + 2b2(0i - 02 - 20) + baa©] 

* [cl(0i - 6) + C2(02 + 0) - (.Ct.. ’ .,.9.2)eo3/^^ 

(1- q)Z 

G2 = -[2bi(0i + 02 ) - 2b2(0i - ©2 - 20) - b30o] 

= [C2(0i - 0) + cl(©2 + 0) + (C| - Co)0 ^3/it^ 

(1 - q)^ 

The parameter H represents a reduced extenslonal stiffness and 61*62 
account for bending-extension coupling. The two definitions of 6j and 62 can 
be shown to be identities by equations (A9) and (AlO) and must be true in order 
for the matrix t to be symmetric. Equations (A9) and (AlO) can also be used to 
calculate all derivatives with respect to q that are required. 

At this point the form of the t matrix is identical to that of reference 1. 

The only difference is that the definition of H and 6^ contain extra terms as a 
result of the member curvature. The transformations necessary to obtain the 
tangent stiffness matrix are based on the geometry and eqilibrium of the member 
under incremental displacement. The final equations are identical to those in 
reference 1 but are developed in a slightly different manner. The forces on 
the member in global coordinates are given by 

F » TVfp (A19) 
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where 



1 

0 

0 

-1 

0 

0 

B - 

0 

1/L* 

1 

0 

-1/L* 

0 


0 

1/L* 

0 

0 

-1/L* 

1 


The Incremental forces due to a small change In the displacement vector are 


given by 


6F » TTBT6f^ + C«D 


where 


C . [’lllBl) f , f f ) 

^ r 3D? r 3Dfi r^ 


Using equation (A15), equation (A21) can be written as 
6F * rVt 6dp + C «D 

The incremental relative displacements are related to the Incremental local 
displacements by 


6dp B6d 


which can be written In terms of global displacements as 


6dp = BT6D 

Performing the differentiation Indicated In equation (A21) and using equation 
A(24), equation (A22) can be written as 

6F « TT(BTtB ♦ 6)T 60 


(A20) 


(A21) 


(A22) 


(A23) 


(A24) 


(A25) 
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where 


and 



symmetric 


0 

-Q 

0 

0 


-Q 0 
P 0 
0 0 
Q 0 
-P 0 


0 


Q = (mj + m2)/L* 


(A26) 


(A27) 


Thus the global tangent stiffness can be obtained from the conventional coordinate 
transformation oerating on a matrix k defined as 


k = B^tB + G 


(A28) 


The matrix k is symmetric with the independent elements given as follows 


kii = ki4 = k44 = EA 
LH 


^12 = ki5 = k45 = -k24 


•^13 = -1^34 = IJGi 
H 


EA (Gi + G 2 ) + 


mi +_mp 


kl6 

k22 

^26 


-k46 = EAGo 


•^55 " -^25 =^C2(ci+C2) U + ^ (Gi+G2)2] - P 

*^35 “ i [(C 1 +C 2 ) n + EAL 6 i(Gi+G2)] 

L* IT 

-k56 = 1 C(ci+C2) El ♦ EAL G2(Gi+G2)] 

"C~ H 


P 

L* 


(A29) 
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^36 “ ^2 £L 

= CiU* ^ G2 
L H 

The sequence of solution steps when 63 and 03 -u^ are known Is to solve 
equations (A7) and (A8) for pj which can be done Iteratively. Then mj and m3 
can be determined for equation (AS) and the stiffness matrix k from equation 
(A24). If the Interior deflection of the member Is desired. It can now be 
determined rather easily from equation (A4) using the following expressions 
for the constants of Integration. 

Ai = -A3 = mj/pi 
^2 = (mj + m3)/pi 
A4 = -A3 + L(6j- B) 

APPENDIX B 

BUCKLING OF LATTICE STRUCTURES WITH IMPERFECT MEMBERS 

For certain structural configurations having the proper symmetry, the 
initial ro'ilinear state may be rather simple consisting of shortening of indl- 
vldiia; members due to mutual restraint of adjoining members no end rotation 
^ccurs. Examples of such structures are shown In figure 7 which was taken from 
refer 3nce 5 where the buckling of these configurations was studied. The cylin- 
drical con F-; durations when compressed axially, shorten without bending even If 
Individual members are not straight provided all similar members have the same 
Imperfection. Of Interest may be the load level at which such a structure 
bifurcates from this Initial state. The tangent stiffness matrix for the 
Imperfect member with Bj = 63 = 0 could be used In a conventional buckling 
analysis to determine this bifurcation point. Considerable simplification for 


(A28) 
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this case Is possible which will be Illustrated by giving the changes to the 
analysis .>f reference 5 which are required to treat imperfect members* 

The analysis of reference 5 applies to a repetitive structure so that any 
Imperfection must be Identical for each group of members. For this analysis 
one Imperfection, a displacement normal to a cylinder containing the nodes, Is 
assumed to exist In the diagonal members and another In the ring members. 

Only the upper half of the 12 x 12 three dimensional stiffness matrix Is 
required In the analysis so the stiffness terms associated with deflections 
are taken In the plane of the Imperfection from the present planar analysis. 

In reference S the upper half of the complete member stiffness matrix Is 
partitioned Into two 6x6 matrices R and S. The following are the changes 
required from the usual representation In order to account for the Imperfection 


Rll « -$ii « EA/LH 

(Bl) 

R55 • CiEl/L + EALG{/H 

(B2) 

$55 « C2EI/L - EALGf/H 

(B3) 

•^15 • -Sl5 * -S 5 I 

(B4) 


If there were no Imperfection H>1 and Gx>0 and equations (Bl) to (B4) 
are Identical to the standard expression for straight members used In reference 
5. The global stiffness matrix K which yields the buckling eigenvalue given 
In reference 5 In terms of the elements of R and S Is still applicable for 
the modified terms given In equations (Bl) to (B3). However, In reference 5, 

R 15 , $ 15 , $51 which are defined by equation (B4), were all zero so that additions 
to the global stiffness matrix K must be made as follows 

aKi 4 • [4$i5SeCeScS„3d 

4Ki5 * [4(Ri5 ♦ Si5CgC^)cJc^3<j 
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AKjg » [-^Si5c5sj,SgS^]<j 

AKj, ■ [(-4(Ri ♦ S,5CjC„)sJc„]j ♦ [-2(R,5 * S,5C,)c„]r 
AK 25 * C-^Si5C0Se4sgS^3jj 
AK 26 • C-4Si5CeSeCj,SeS„](l 
AK34 » [^SisSgS^CgS^lj, 

AK35 ” C“^SjgC0S0Cg|S(|S£C,^3(l 
AK36 * C4Si5C0S0sJsgS^]jj 

where aK^j Is the addition that roust be roade to K^j with the notation the saroe 
as reference 5, 
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Table 1. The Elastica Problem of Figure 4. 


f/pg , 

ea 



y/L^ 

y/L^ , 

y/b^ 












degrees 


1.015 

.970 

.971 

.970 

.220 

.212 

.218 

20 

19.37 

19.86 

1.063 

.881 

.888 

.881 

.422 

.110 

.422 

40 

38.90 

39.94 

1.152 

.741 

.757 

.741 

.593 

.576 

.593 

60 

58.35 

59.98 

1.213 

.560 

.589 

.558 

.719 

.699 

.719 

80 

n.7i 

79.99 

1.518 

.349 

.398 

.348 

.792 

.770 

.790 

100 

96.97 

99.98 

1.884 

.123 

.198 

.125 

.803 

.788 

.800 

120 

116.12 

119.95 

2.541 

-.107 

-.001 

-.102 

.750 

.753 

.745 

140 

135.14 

139.86 

4.029 

-.340 

-.208 

-.324 

.625 


.617 

160 

154.05 

159.69 

9.116 

-.577 , 

-.459 

-.529 

.421 


.426 

176 

169.47 

175.63 


1 Reference 3 

2 No Interior node 

3 1 interior node 


20 















DOUBLE LACED CYLINDRICAL TRUSS 


ORIGirML PAGE fS 
OF POOR QUALITY 





O SPACE FRAME ANALYSIS 


ORIGINAL PAGE IS 
OF POOR QUALITY 





1 


